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In the framework of the QCD dipole model at high energy, we present an analytic evaluation 
of the dipole pair density in two limits in which the parent dipole is much larger /smaller than 
[ — . the distance between the two child dipoles. Due to conformal symmetry, the two limits give an 

identical result. The power-law correlation between dipoles explicitly breaks the factorization 
of target-averaged scattering amplitudes. 
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1 Introduction 



A hadron in the infinite momentum frame is a complicated system of small-x gluons. 
While the energy evolution of the average gluon number can be described by the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) equation [1], the wavefunction of a hadron contains more 
information than just the average number. For example, the fluctuation of the gluon 
number plays a crucial role in the evolution of scattering amplitudes towards the unitarity 
limit, and has recently attracted considerable interest [2-9] in the context of saturation 
physics [10-12]. 

Another important characteristic of the hadron wavefunction is the correlation of gluons 
in the impact parameter space. In the dilute, non-saturated regime, soft gluons are nec- 
essarily correlated because they originate from a common ancestor via gluon splitting. 
The process can be most easily described in the QCD dipole model formulated in the 
large N c approximation [14-16]. In this approach, the evolution of the 'parent' dipole 
(a quark-antiquark pair) proceeds via dipole splitting with certain probability computed 
in perturbation theory (Fig. 1). Since the probability depends nontrivially on transverse 
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Fig. 1. A cascade of dipole splitting. Thin lines represent child dipoles produced by the parent 
dipole (xq, x\). 



coordinates, 'child' dipoles will be distributed in the transverse plane with characteristic 
correlations between them. Although this dynamics is built-in in the numerical Monte- 
Carlo simulation of this model [2, 17], so far there have been only few analytical insights 
[15, 16, 18, 19]. (See, also, [20].) In this paper we evaluate the dipole pair density in cer- 
tain limits and find the power-law correlation between dipoles at large distances with 
the power determined by the conformal weights of the BFKL eigenfunction [21]. As an 
immediate consequence of our result, we shall show in Eqs. (3.18) and (4.21) that the fac- 
torization of dipole scattering amplitudes is explicitly violated by a position-dependent 
multiplicative factor 



{T{T 2 ) wci 2 (Ti)(T 2 ) c 12 >l, 



(1.1) 



where T is the single dipole scattering amplitude and (...) denotes the averaging over 
the target wavefunction. Eq. (1.1) should be contrasted with the fact that scattering 
amplitudes computed in the BK-JIMWLK framework [13] essentially factorize 



(7\T 2 ) » (TiXT; 




(1.2) 



The gluon splitting diagrams which lead to Eq. (1.1) are not included in the BK-JIMWLK 
equation which rather sums gluon recombination diagrams. Thus it is not surprising that 
the factorization in Eq. (1.2) does not hold for a more general evolution. While Eq. (1.2) 
may be valid if one starts with a large nucleus with totally uncorrelated partons [12] 
and follows the BK-JIMWLK evolution up to not so high energies, it is likely that the 
correlation in the transverse plane developed in the dilute regime significantly affects the 
nonlinear evolution of hadrons as in the case of the gluon number fluctuation [3-8] . 
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2 Single dipole density 



In this section, we review the properties of the single dipole distribution. The techniques 
used here can be directly applied to the analysis of the dipole pair density in the next 
section. The single dipole density evolved up to rapidity Y is given by 

n Y ( X01 ,x) = 2 X 4 f^e^ Y ^, (2.1) 
x z J 2m 

x oi = x — xx and x = x 2 3 = x 2 — x 3 denote the coordinates of the parent dipole and the 
child dipole, respectively. We shall use the letter x for two-dimensional real vectors and 
z for corresponding complex coordinates. By slight abuse of notation, we use x also for 
the magnitude of two dimensional vectors, x is the usual BFKL eigenvalue 

X (7) = 2a s Re{^(l)-^(7)}- (2.2) 
The saddle point of the 7-integral is given by 

X '( 7 )F = ln4 (2-3) 
x A 



When #01 > x, the saddle point is in the region 1 > 7 > |, and 

ny(x i,x)~e*W y ^y \ (2.4) 
n is proportional to the scattering amplitude between dipoles of sizes x i an d x. 

TXO?X^ 

T Y (x i,x) = s —n Y (x i, x) (2.5) 

27 2 (i -7) 

Eq. (2.1) is integrated over the impact parameter b = X2 + X3 — ga±gi between the parent 
and child dipoles . The ^-dependent distribution is 

x Jf uE ™-\b+^-u,,b-^-u,)E^(?f- U ,-?f-<j). (2.6) 
E is the eigenfunction of the SL(2,C) group 



h / _ \ h 
ZOI \ 



^ ^07^17 / 

E h ~ h *(x 0y , x l7 ) = E 1 -*' 1 -*^ x ly ) (2.7) 
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with h = + iv, h = ±±2 + iv = 1 - h*. Eq. (2.1) is obtained from Eq. (2.6) by 
integrating over b 

n Y (x m ,x) = J d 2 b n y (x m ,x,b), (2.8) 

keeping only the n = term and identifying h — \ + iv = 7. The w-integral in Eq. (2.6) 
has been carried out in [22,23]. Due to global conformal symmetry, the result depends 
only on the anharmonic ratio 

P - (2.9) 

Z02Z13 

The n = term gives, 
2z/ 2 / 

no^xoi, x, b) = — VIP| 2(1_7) ^(1 - 7, 1 - 7, 2(1 - 7); - 7, 1 - 7, 2(1 - 7); p) 

+&* ) Jp|^F(7,7,27;p)F(7,7,2 7 ;p)), 

(2.10) 

where F is the hypergeometric function and 

2 4 - r(|-^)r(i + ^) 
= * 17 r(§ + «/)r(i-«/)- (2 - n) 

Consider the case x i ^ x and look at the region of small impact parameters xqi ^ 111 
this region, 

— 4z 

, |p|<l, (2.12) 

and one may approximate F(...; p) ~ 1. We obtain 1 

ny(xoi,x,6) = / pn Axonx,b)e*W « i~ / ^ ^^^ f^V " + c ' c - 

(2.14) 

Comparing Eq. (2.1) and Eq. (2.14), one sees that in the saddle point approximation, 

n Y (x 01 , x) ~ x 2 0l n Y (x i,x, b < x i). (2.15) 



1 Ref. [24] uses the following approximation 

16 2iu 

b ,v^K 3 ^. (2.13) 
iv 

/16a; 2 

This is valid as long as v is close to zero and leads to a factor ( J? 1 J . However, in our 
case the saddle point for v is not assumed to be small, but rather determined from external 
parameters (dipole sizes). 
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Therefore, roughly child dipoles are uniformly distributed inside the area (c.f., Eq. (2.8)). 
On the other hand, the dipole density at large impact parameters b ^> x 0i are suppressed. 
Indeed, in this region, p pa and 

1 r u 2 h* St 2 t 2 \ 7 
n Y (x 01 ,x, b)^^ 2 j d/-^- pg^J e^ Y + c.c. (2.16) 

At the saddle point, n(b) ~ 1/6 47 where 7 is determined from x'(l)Y = In Jt x2 ■ 

Let us compare this 6-dependence with that of the saturation momentum. The dipole- 
dipole scattering amplitude at a fixed impact parameter b is 

T Y (x 01 ,x,b) = J d%'-^A (x,x',b-b')n Y (x 01 ,x',b'), (2.17) 

where A is the dipole-dipole scattering amplitude in the two-gluon exchange approxi- 
mation. Since A (b — b') decays like 1/(6 — b') 4 , one may approximate 



/f d x 
d 2 b'A (x,x',b-b') J - — -n Y (x 01 ,x',b) 

= ^a 2 s J _|_ r |(i + l n !>) ny ( Xoi;X ' )& ) (2.18) 



where r< = min{x, x'} and r> = max{a;,2;'}. Using the large b form of n, Eq. (2.16), one 
obtains 2 



na 2 x 2 



T Y (x Q1 ,x,b) pa -^—^ n Y {x Q1 ,x,b). (2.20) 



The local (6-dependent) saturation momentum can be determined by the constancy of 
the exponential factor of T in the integral representation along the line x = l/Q s (b,Y) 
[25, 26], and reads 



Ql(bX)~ X f^ Y , (2.21) 



where 7 S pa 0.628 solves x'ils) — Repeating the same procedure for x m 3> b, we get 

Q 2 s (b,Y)r,^ e ^ Y . (2.22) 
x 01 



2 Compare with Eq. (2.5). The factor 2 difference in the denominator is due to the definition 

T Y (t = 0) 



= J d 2 bT Y (b). (2.19) 
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Fig. 2. Graphical representation of Eq. (3.1). 



If we take x to be close to the saturation line ~ l/Q s (b,Y), 7 ~ 7 S , and the geometric 
scaling [25-27] holds locally in the two (6>i i, b <C x 01 ) regimes [28] 



Ty(x i,x,6) ~ x n Y (x 01 ,x,b) ~ (x Q s (b, Y)) % 



(2.23) 



3 Dipole pair density 



In this and the next section, we analyze the dipole pair density rS 2 ^ in two different ways. 
We start with the exact expression for the pair density as derived in [18] (see, also, [29]). 



n { Y ) (x 01 ,x aoai ,x bobl ) = / dhdh a dh b 



1 



2x 2 ,^„,x\^ Jo 



Y 



dy e 



x (h)y+(x(ha)+x(h b ))(Y-y) 



-'a ai J 'bof>l 

x / d 2 x a d 2 x l3 d 2 x 1 E h ' h (x 0l , x ll )E h ^ ha (x aoa , x aia )E hb ' hb (x bo /3, x blP ) 



X 



I d 2 X2d 2 x 3 d 2 ^ e hm {x2 ^ X3i )E h ^*(x 2a , x, a )E h ^{x^ X4fi ), (3.1) 

J Xo^XqaXa') 



where (x ao ,x ai ) and (x bo ,x bl ) are coordinates of the child dipoles of interest. We intro- 
duced a compact notation 



dv 



2u 2 + n 2 /2 



7P 



(3.2) 



A graphical representation of the coordinate integrals is shown in Fig. 2. We will be 
interested in configurations where the two child dipoles are small (typically of the order 
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(A) (B) 

Fig. 3. Case (A): All dipoles are small compared with their mutual separations. Case (B): Child 
dipoles are deep inside the large parent dipole. 



of the inverse saturation scale ~ 1/Q S ) an d far away from each other, 

_ X ao ~r %ai x b x bi „n 

%ab %a x b = ^ ^ 2-ooai, %bobn I" - "/ 

and try to extract the leading x a b dependence of n^ 2 \ This leaves us with two interesting 
(and in fact, tractable) situations (see, Fig. 3): (A) The parent dipole is also small Xqi ~ 
^0001,^6061 < x ab- (B) The parent dipole is large x 01 > x ab > x aoai ,x bobl . 

Let us first consider the case A. In the next section we will discuss both cases in a unified 
way. The last line in Eq. (3.1) is the triple pomeron vertex in perturbative QCD [30] at 
large N c . It has the form 

1 

f(h, h a ,h b ) 



(z a - z l3 ) 1+h - h o- h b(z{ 3 - z 1 ) 1+h ^- h b- h (z 1 - Za y+h b -h-h a 

1 

(Z a - Zp) l + h -~ h "-~ h b(zp - Z^Y+^-^-h^ - Z a y+ h b-h-h a ' 

This structure follows immediately by noting that the last line of Eq. (3.1) and Eq. (3.4) 
transform in the same way under the SL(2,C) transformations of z a , zp, z 1 

z ^az±j3 ( a5 _p 1 = l y (3.5) 



The function f(h, h a , h b ) can be found in [31, 32]. Next we turn to the remaining integrals 
ePxaePxpePx^ in (3.1). Since all dipoles (parent, children) are assumed to be very small, 



typically x 0a , x b, x a b ^> x i, x aoai , and we may make approximations 
[ Zm \ ~ [ Zoi \ ( Za ° ai \ ~ [ Zat > ai \ ( Zht)bl \ ~ ( Zb » bl \ 

\ z 0j z l-y) \ Z 0-yJ \ z a a z a ia J \ z l a J \ z b a l3 z b x f3 J \ z b0 J 



We will see later that with this replacement one makes a mistake in the overall factor of 
by 8. After this approximation, we are left with the integral 



/ 



j2„ .2„ 1 



d x a d xpd x 



7 ~2h„2h a -y 2hb v 1 +h-h a -h b l+h a -h b -h l+h b -h-h a 
Z 0~f Z aa Z b(i Z aj3 Z (3-y Z 7" 



* ~2h ~2h a v 2h b -l+h-ha-h b -l+h a -h b -h-h b -h-h a ' (^-7) 
Z 07 z aa Z bf3 Z a(3 Z /3~/ Z ~i a 



One can check that this integral transforms in the same way under the SL(2,C) transfor- 
mation of zq, z a , Zb as 

(3.8) 



h+h a -h b h+h b ~h a h a +h b -h h+h a -h b -h+h b -h a -h a +h b -h 
-Oa ^06 z ab z 0a Z 0b Z ab 



The coefficient can be easily obtained. In the dominant case of n = n a = rib = where 
h = h = | + iu = 7, h a = h a = 7a , h b = hb = 7 &, (Generalization to the case h ^ h is 
straightforward.) 



J d 2 x a d 2 xpd 



= g(l,la,lb) 

" | z Qa\ 2[l+la " 7b) | z 0b 1 2 ( 7 +Tf-^) | Z ab \ 2 ^ + ^~^ 

(3.9) 

When 7 > 1/2, a pole at z 7 = zq is not integrable. The following result should be regarded 
as analytic continuation from convergent values of 7's. Using a conformal transformation, 
one can set z a — 0, Zb — 1, zq — 00. 

r 1 1 

<?( 7 , 7 a,7b) = J d 2 x a d 2 xpd 2 x 1 ^ 4 ^ i _ | Za/3 |2( 7 - 7a - 76+ i)| Z/g7 | 2 ( 7a - 7b - 7+ i)| %Q |2( 7 ^ 7 - 7a+ i) 

(3.10) 

Evaluating the integrals in the order of d 2 x 11 d 2 xp and d 2 x a , one obtains 

^(7, 7a, 7b) = 7T 



3 r(i - 2 7 )r(i - 2 7o )r(i - 2 7fe )r( 7 + 7a + lb - 1) 



r(2 7 )r(2 7a )r(2 76 )r(2 - 7 - 7a - 76 ) 
x r( 7 + 7a - 7& )r( 7 + 7fc - 7a )r( 7a + 7fc - 7 ) n 
r(i + 7fe - 7a - 7 )r(i + 7a - 7b - 7 )r(i + 7 - 7a - 7 b) ' 
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and therefore, 

= f didr( a dr( h 1 2 3(7, 7 Q , 7&)/(7, %, 7b) / ^ 
/ d-fd-f a d-f b — — p(7,7a,76)/(7,7a,76) / 



n = / ^7^7. = 2 2 g(7, 7a, 7b)/(7, %, 7b) / <*y 2(^.-^) 2 (7+7 »- 7n ) 2( 7a+7t - 7 ) 

aoai^booi J-Oa x Qb x ab 



J aoa\ bobi 

x exp ( X ( 7 )y + ( X ( 7a ) + x(7b))(y - *) - 7 In ( ^) - 7. In (j^-) - 7o ln ( 

(3.12) 



The remaining integrals may be evaluated in the saddle point approximation. The saddle 
points for y, 7, 7^ are given by the solution to 

X(l) = x(la) + X(lb), (3.13) 

„2 ^2 



X'(7a)^-y) = ln-5^»l, (3.14) 

x 0b x a a 1 



xXlb){Y-y)=\n^^^l, (3.15) 
X '( 7 )l/ = ln44» 1 J (3-16) 

X ab x W 

and the dipole pair density behaves like 

e {xi~la)+x{~lb))Y xllx 2la xVl 

n Y {x i,x aoai ,x bobl ) 2( 7+7o - 76 ) 2( 7+T6 - 7a ) 2( 7o+76 - 7 )- 

a oai &obi x 06 x ab 

The Xqi behavior of was pointed out in [16]. (See, Eq. (A. 2) of [16].) From Eq. (3.16) 
we see that \ < 7 < 1, and this justifies the conjecture below Eq. (A. 7) of [16]. The factor 
l/ x la (with 7 a = 76 = \) was found in [19] in the context of dipole production at large 
transverse distances. 

The factor x^ la+lh ^ characterizes the correlation of dipoles in impact parameter space. 
To see the significance of this factor, consider scattering of two dipoles x aoai ,x bobl on a 
target dipole X01 at large impact parameter. The scattering amplitude is given by 

d 2 x d 2 x' 



T Y (x i;x aoai ,x bobl ) = J 77^77^^ / d 2 bd 2 b'n ( y\x 01 ,xb,x'b') 



At X + Xi , , X0 + X1 , 

x A(x aoai x, x a b)A(x bobl x , x b b ) 

^afyloa^bob! (2)/ x 

, " V v' 01 • •' (':,'/; • •' M'l ) 



16 7a 2 (l-7a) 2 7 6 2 (l-7 b ) 2 ' 

f x 0a x 0b\ 2 ( 7a+7b_7 ) 

TY(x Q1 ,x aoai ,x 0a )TY(x 01 ,x bobl ,XQ b ) — - — . (3.18) 

\ x 0l x abS 
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In the second line, we have used the same approximation as in Eq. (2.18). The last line 
should be taken with care since 7 a as determined from Eq. (3.14) does not coincide with 
the anomalous dimension of T(rr i, x aoai ,Xo a ), the latter being determined from x'(la)Y — 

4 

In a " a (c.f., Eqs. (3.14), (3.16) and note that 7 ^ j a ). Even if we neglect this difference, 

X a a 1 x 01 

we see from Eq. (3.18) that the factorization of two-dipole amplitude is explicitly violated 
by a nontrivial position-dependent factor. 3 

\X lX ab J 



4 Improved calculation 



Let us return to the integral appearing in Eq. (3.1). 

I = J d 2 x a d 2 xpd 2 x 1 E h ^ h {x 0l) x ll )E h ^ ha {x aoa) x aia )E hb ' hb {x bo p,x bl ( 3 ) 
x / d2x f^ X4 E h ^(x 2j , x 3l )E h ^*(x 2a , x 4a )E h ^*(x 3f} , X4P ). (4.1) 

J XoiXiaXa') 



rt* 1 -' rv*'-' ry*- 1 

X 23 X 34 X 42 

Instead of first integrating over £2,3,4 ('reggeon coordinates') as we did before, now we 
integrate over x a ^ a ('Pomeron coordinates', see, Fig. 2) first. 

'=7£3»/ ^ X /j€ Xl (WWl, 7. 27; K)F(i, 7, 2 7 ; p) + c.c 

[ZTl ) J x 23 x 34 x 42 ^ 

X (b^uJPal^Fi-fa, 7a, 27a! Pa)F(-f a , ~f a , 2^ a ] p a ) + C.c)j 

x (b* 0tUb \ Pb \ 2jb F( lb , 7b , 2 7b ; p b )F{ lb , 7b , 2 7b ; p b ) + c.c.) , (4.2) 

where 

^01^23 Z aoai Z24: Z bobl Z 34 /. n 

P = , Pa = , Pb = , (4.3) 

^02^13 Z ao2 Z ai4 Z bo3 Z bl4 

are anharmonic ratios. By assumption, p a and p b are small, and one may approximate 
F {■■■', p a ,h) ~ 1- A remarkable point is that p is small both in the limits of Xoi ~~ 00 an d 
x i — > and one may approximate F(...;p) m 1. Expanding the brackets, we get eight 
terms. The first term reads 



& * b * f d 2 x 2 d 2 x 3 d 2 x 4 / a:oiX23 \ 27 / x a(m x 24 \ 2 ^ / x bobl x 34 \ 2 ^» 

7T 6 °' V °'" a 0,1/6 J X 2 23 X 2 34 X 2 42 \X 02 X 13 ) V X a2 X a4 J V X b3 X b4 J 



3 We note that the two-dipole scattering amplitude that appears on the right hand side of the 
BK-JIMWLK equation [13] is for contiguous dipoles, x ai = x bo . Our present approach does not 
apply to this interesting case since we assumed x ab 3> x aoai , x bobl . 
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If we take the limit x 01 — > 0, this is the same integral which gives the triple Pomeron 
vertex /(7,70,7ft)- 



/• d 2 x 2 d 2 x 3 d 2 x 4 fxm^\ 2 "* / x aoai x 2 4 \ 2la / ^6061^34 \ 276 

7 X2 3 X 34 X 42 VX02X13/ V aJa2^a4 / V ^63^64 / 

f d 2 x 2 d 2 x 3 d 2 x 4 / x 23 \ 2j ( x 24 \ 2la ( ^34 \ 

W ^23^34^42 VX 2^03/ \X a 2Xai) \X b3 X b4 J 



2~f b 



27„27 a 27(, 
~ ^-Ol^aoai^feofc 

2 7 2 7a 27 6 /(7; 7a; 7b) /- r\ 

" ^Ol^aoai^ftofcl 2(7+ 7o - 76 ) 2(7+ 76 - 7a ) 2( To +7 6 - 7 ) ' ^ 4 -°J 
x Oa x Ofe X afe 

For a later use, we note that when r ) a = lb, 

h = >£,/h,«J {^-f (^f , (, 01 0). (4.6) 
Eq. (4.5) should coincide with our previous result Eq. (3.12), so we obtain an identity 

^(7,7a,7fc)/(7,7a,7b) = ^ b 0,u b *0,u a b 0,uj(l,la,lb)- (4.7) 

Eq. (4.7) is a straightforward generalization of the relation between /(7, 7 a ,7&) and 
/(7,7 a , 7b) derived in [32]. We also see that the previous approximation Eq. (3.6) misses 
the seven other terms in Eq. (4.2) which contribute equally to due to the symmetry 
7 <-> 1 — 7 of the integrals / d"fd^ a d%. We take this into account by multiplying Eq. (4.4) 
by 8. 

T ]_,* ,* ,* f d 2 x 2 d 2 x 3 d 2 x 4 f x 01 x 23 \ 2 ^ f x aoai x 24 \ 2 ' ya f x bobl x 34 \ 2 ^ 

A W «J x| 3 x| 4 xi 2 U2W UmsmJ 1 J 

Now we would like to evaluate this for xoi — > oo. In the following, we assume that 7 a = 7&, 
which will be approximately valid at the saddle point when x aoai ~ x bobl . Then Eq. (4.8) 
takes the form 

1 /™ ™ \ , j2„ j2„ j2„ /„ „ x <2-y / „ „ „2 \ 2 7a 

!,* !,*2 / x aoai x 6obi 



7 ~ t/oA^ r 
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/• d 2 x 2 d 2 x 3 d 2 x 4 / xoiX 23 \ 27 / 3^24^34^ \ 7 ° „s 

7 ^23^34^42 ^02X13/ \X a2 X a4 X b3 X b4 J 



It is easy to see from the SL(2,C) invariance that the result of the integration must have 
the structure [In fact, this property holds only for 7„ = 7^.] 



I = -S, b t a i X -^^T h (p',p') , (4.10) 

71 \ X ab ) 



where 



p' = (4.11) 

ZQ a Z\ b 
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is the anharmonic ratio of the external points. In the limit x m — > 0, p' — > and /i should 
reproduce Eq. (4.6). 

%',p)^/(7,7a,7a)|pf 7 , (p'^0). (4.12) 
Our observation is that the limit xm — > oo also leads to p' — > 0. Therefore, when xoi — > oo, 

= 8^(7,7a,7a)/(7,7a,7a) fes*^) ^ f^l) ^ , (4 .1 3 ) 

and we obtain the behavior of at the saddle point 

1 / \ 2 T a / \ 27 

( 2 )/' „. .„ „ „ „ \ 1 ^xila)Yi* 7*2 /•/ „, „, \ / ^aoai^feobl \ / x 01 x ab \ 



^aoai^feobl V X ab / V^Oa^lft/ 

(4.14) 



with 7 and 7 a (and also y) determined from the saddle point equations 

Xil) = Mia), (4.15) 



2 „2 



X'(7)y = ln^»l, (4.16) 



^y»2 ,yj 2 

x 01 x aft 

X / (7a)(i r -y)=ln-^-»l. (4.17) 



Note that | < 7 < 7 < 1. 

Eq. (4.14) is valid both for x i <C £ a & and x i ^> rr ab . Due to conformal symmetry, the two 
cases of Fig. 3 are mathematically identical. In the latter case, if the dipole x ab is deeply 
inside the parent dipole x i, one can make an approximation 



\XoaXibJ \XoiJ 

It is interesting to note that in the same limit we may rewrite as 

i-Y 



(2) _ /(7,7a, 7a) 



p 1 

X X X 

\ dy 2~ 7ly (Xoi , iCab) 2 ^Y—y\ x ab, x a^a\ ) 2 ^Y—y[ x abi x bobi ) • 
JO X m X„ h X nh 

(4.19) 



2 2 

^aoai^feobl 1/0 ^Ol ^ab ^ab 



Eq. (4.19) provides an intuitive understanding of the result. The parent dipole x 01 emits 
a child dipole x a b inside the area with uniform probability (c.f., Eq. (2.15)). The 
geometrical factor x 2 ah jx\ x specifies the location of the dipole x a t,. Then the dipole x a b 
splits into two dipoles of similar size ~ x a b through the triple Pomeron vertex, 7(7, 7 a , 7 a ). 
Finally, each of the two dipoles emits a child dipole of size x aoai (or inside the area 
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~ x\ b again with uniform probability, and the two child dipoles x aoai and x bobl roughly 
fall within a distance x a b (see, Fig. 3B). Another representation of is (c.f., Eq. (2.5)) 

r$ oc / dyTy(x 01 ,x ab )TY-y(x ab ,x aoai )T Y -y(x ab ,x bobl ), (4.20) 
Jo 

which may be a useful form to include effects beyond the BFKL evolution. 

Finally, we consider how the approach to saturation is modified in the presence of power- 
law correlations in the target. The scattering amplitude of two dipoles off a large onium 
of size xoi at small impact parameter can be computed similarly as before 

Ty\x 01 ; x aoai , x bobl ) ~ o^x 2 aoai xl obi nY{xQi]x aoai ,x babl ) 

/0; 01 \ 2 ( 2 7a-7) 

~ T Y (x 01 ,x aoai ,b « 0)T Y (x 01 ,x b{)bl ,b w 0) I — ) . (4.21) 

The power law correlation Eq. (4.21) is remarkable in view of the fact that the single 
dipole distribution has essentially no 6-dependence deep inside the dipole xoi- Due to the 

enhancement factor, (j^J ° ^ 1, the problem of unitarity for is severer than 
that for T. The condition T( 2 ) < 1 is roughly equivalent to requiring that the exponential 
factor of vanishes along the saturation line x aoai = x bobl = 1/ Q pa irsat 

2 

2 X { la )Y - 7 In X -f - 2 la HxlQlairsat) = 0. (4.22) 

X ab 

Solving this equation with the conditions Eqs. (4.15)-(4.17), one finds that 

= ^ e 7a # = ^ exp — ir Y ■ ( 4 - 23 ) 

X 01 \ X abJ X ab VTM" 7 / 

Since 7 and 7 a depend on coordinates, there is not a unique way of writing Qpairsat- The 
first expression shows that, when 7 a pa j s , the onset of unitarity corrections is much 
earlier than the single dipole scattering case, Q pa irsat 3> Qs (c.f., Eq. (2.22)). This is 
quite a contrast to the result Q s ~ Qpairsat which would follow from the assumption of 
factorization pa T 2 . The second expression of Eq. (4.23) emphasizes 4 that Qpairsat is 
not an intrinsic quantity of the target, but depends rather sensitively on the configuration 
of the projectile. 
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